Characterizing order in amorphous systems 
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We measure and compare three correlation lengths proposed to describe the extent of structural 
order in amorphous systems. In particular, the recently proposed "patch correlation length" is mea- 
sured as a function of temperature and fragility and shown to be comparable with other measures. 
In addition, we demonstrate that the patch method also allows us to characterize the symmetries 
of the local order without any a priori knowledge of it. 



Amorphous solids are ubiquitous, and include glasses, 
sand piles, emulsions, and colloidal systems [THS]- Such 
systems appear at best to possess short-range order [3J 2], 
in stark contrast to crystals and quasicrystals, whose 
order may well extend over an entire macroscopic sam- 
ple. One striking common feature of amorphous systems 
is slow dynamics: thermal systems become increasingly 
sluggish upon cooling [1 , 3 , showing a remarkable rapid 
rise in viscosity, while athermal systems jam when com- 
pressed [H 13 |S] . While it is easy to understand how mo- 
tion in an ordered system is difficult, it is harder to ratio- 
nalize such slow dynamics in amorphous systems, which 
leads to a fundamental question: Is the extraordinary 
slowing-down related to some subtle growing structural 
order? This question is the subject of much current activ- 
ity [3J HI ITrfT^] , and its resolution would lead to a better 
understanding of glassy phenomena, and could facilitate 
the testing of various theories of the glass transition [3J. 

Classical methods used to probe structural order show 
no dramatic change in the structure of a material as it 
approaches the glass transition [3J. These methods are 
based on scattering experiments and only give access to 
two-point correlations, such as the pair correlation func- 
tion and the structure factor. Thus, if there is any ("hid- 
den" ) order present in glassy systems, it would have to be 
very subtle, and inaccessible to conventional techniques. 

A major obstacle in measuring the extent of order in 
conventionally amorphous systems is our lack of knowl- 
edge of what that order might be. This is unlike more 
familiar systems in condensed matter where we identify 
the order parameter at the outset. Can we, in the ab- 
sence of a well-defined order parameter, define a charac- 
teristic correlation length? Since the nature of a possible 
order in a glassy phase is not known, a good definition 
of a correlation length should assume nothing about the 
nature of the order it seeks to quantify. That is, struc- 
tural correlation lengths should be defined in an order 
agnostic way. Such a procedure would equally well aid in 
the identification of possible hidden order in other, non 
glassy systems, such as colloidal suspensions. 

Interest in the glass transition has spurred several re- 
cent proposals of structural correlation lengths in super- 
cooled liquids [3J [TTJ [T2]. For example, the "point-to- 



set" length [H [3 [8] assumes nothing about the nature 
of the developing order, but has the disadvantage of be- 
ing difficult to implement experimentally since it requires 
freezing a portion of the system in its equilibrium con- 
figuration while letting the rest evolve. The proposal of 
Ref. 1 11 j relies on a-priori knowledge of the developing 
order, which is the case for the frustrated system studied 
here. Although it is not order- agnostic, we use it in this 
Letter as a benchmark for comparison. 

Recently, a novel, generic method to define and mea- 
sure a static correlation length [9l [10] has been pro- 
posed, based on the frequency of occurrence of motifs 
("patches") of various sizes. The main idea is to com- 
pute the entropy of patches of a given size R appearing 
in a system [15] . and then to look at how this entropy 
scales with R. This scaling yields information about the 
extent of order in the whole system, even if no order pa- 
rameter has been identified, and allows us to define a 
"patch correlation length" . 

There are two main points to this paper: First, we cal- 
culate the patch length for the frustrated atomic system 
of Ref. [llj as a function of temperature, and compare 
it to other proposed correlation measures. Second, we 
demonstrate how patch entropy considerations may be 
used to determine the nature of the developing order. 

Although the patch entropy method does not require 
patches of any specific shape, we shall find it convenient 
in what follows to consider spherical patches of radius R. 
To get a feeling for the behavior of the patch entropy, 
we first consider a perfect undefected crystal. It is easy 
to check that for R greater than some radius (roughly 
the size of the unit cell) , the number of possible patches 
is constant, independent of R, and therefore so is the 
patch entropy £>(-??) [HI HH]- The opposite extreme is a 
totally random system, where the number of patches is 
exponential in the volume V(R) of the patch, leading to 
S(R) oc V(R) - that is, it is extensive in patch size. For 
a quasicrystal, S(R) cx logi?, i.e., the entropy is subex- 
tensive in patch size, albeit greater than for a crystal. 

Next, consider a system ordered only up to a finite dis- 
tance £. It might be expected that in such a case, the 
system would, on short scales, organize into low energy 
patches which are present in the ground state - this is 
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what happens, for example, in a polycrystal. In this case, 
small patches (i? < £) will occur often, and the patch en- 
tropy will be subextensive in R for R < £. On the other 
hand, since any given large patch will occur exponentially 
infrequently, S(R) ~ V(R) for R > £. (For the poly- 
crystalline case, for example, patch congruence would 
require the identical configuration of grain boundaries.) 
The scale £ where the behavior of S(R) changes behavior 
from subextensive to extensive was proposed [HI HO] as a 
measure of the correlation length of the system; this is 
the patch correlation length In practice, it is easier to 
measure this length by measuring the inverse slope of the 
extensive regime of the entropy; this is what we do in the 
following. In particular, one defines a correlation volume 
V* ~ £_g (d being the space dimension) by noting that in 
the extensive regime, S(R) cx V ^ . (Since this is only a 
proportionality relation, we cannot fix the absolute value 
ofCs'.) 

In [10] , it was demonstrated that this proposal for 
defining order applies not only to systems whose ground 
states are crystalline or quasicrystalline, but also to other 
kinds of exotic ground states whose order is not detected 
by scattering. Thus, this method may reveal "hidden" 
order present in a system. As argued in Ref. [10], for 
a non-stationary system evolving in time, the relevant 
patches to be studied are derived from the density pro- 
file which persists after high frequency vibrations are av- 
eraged away. This consideration is based on a separa- 
tion of time scales, which is seen in the relaxation of 
glassy systems: rapid motions occur on the so-called f3 
relaxation time Tp, whereas the long time a-relaxation 
(t q ) corresponds to collective structural relaxation. The 
rapid motions then average out over times r such that 
T/3 <C t <C r a , and we are left with a density profile on 
which we may measure the patch entropy. 

In what follows, we will analyze static configurations of 
the model described in [TTJ that are obtained by quenches 
to low temperatures. The use of static configurations is 
motivated by experimental systems for which snapshots 
are accessible while time averaging is difficult. In this 
case, however, no two patches are exactly identical, and 
we need to define a procedure for quantifying the "dis- 
tance" between two patches. For the purpose of comput- 
ing the patch entropy, we shall consider as congruent two 
patches whose distance is less than some cutoff. Our ra- 
tionale is that two very similar patches from a static con- 
figuration would, could we observe them in time, evolve 
to the same average density profile. In consequence, a sin- 
gle (sufficiently large) snapshot of a system should serve 
in lieu of time averaging. We suggest that choosing the 
cutoff from the value of the plateau seen between the (3 
and a regimes in the mean square displacement in glassy 
systems [Tj should be analogous to averaging over a time 
corresponding to the end of this plateau. 

The detailed procedure is as follows: (i) Select a par- 
ticle i in the system and consider the surrounding con- 



figuration contained in a sphere of radius R centered on 
i - denote this patch IT(i?); (ii) Superimpose this patch 
on ILj(R), the corresponding patch of another particle 
j, so that the centers coincide; (iii) For each particle in 
TLi(R), compute the distance d from the nearest particle 
of the patch Uj(R). If d is less than a given threshold d* , 
consider these two particles to be matches; (iv) Rotate 
the patch by some small angle 6a and repeat step iii. If, 
for some angle a, each particle inside the (rotated) patch 
ITj(_R) matches one in Hj(R), consider the two patches 
to be congruent |16j : (v) Build a congruence matrix C 
whose element c^- equals 1 if IT and IT, are congruent, 
otherwise set — 0. Next, symmetrize C: if cy ^ cji, 
set — Cij = (so that if IT is congruent to IT, then 
the reverse is true as well); (vi) Construct equivalence 
classes by grouping patches such that all the patches in 
a class are congruent with at least one other member of 
the class [17]. Intuitively, patches belonging to the same 
equivalence class look similar, while those belonging to 
different classes look different. 

Having constructed the equivalence classes, we com- 
pute the associated patch entropy S(R): 

S(R) = -Y / P a (R)^P a (R) (1) 

a 

where the sum runs over equivalence classes, and P a (R) 
is the probability of a patch to be in the a th class. We 
next examine the functional dependence of S on R. As 
discussed earlier, one can extract a correlation length ei- 
ther by identifying the crossover length where S(R) goes 
from subextensive to extensive in R, or by measuring the 
slope of the extensive regime. 

We will now apply these ideas to the model glassy sys- 
tem of Refs. [TTJ [T31 [T3], which consists of monodisperse 
Lennard- Jones particles placed on a two-dimensional 
negatively curved surface. The system tends to have a 
local hexagonal order, whose spatial extension is frus- 
trated due to the curvature of the surface. Nevertheless, 
hexagonal or hexatic order persists over a temperature- 
dependent length scale, and may be quantified with an 
appropriate correlation function. In this model, the dy- 
namic fragility of the system can be tuned by varying the 
curvature of the surface [T5] . 

To evaluate the patch entropy of this system, we must 
first determine the relevant value of the threshold d* . As 
argued above, a reasonable choice for d* is the value of the 
mean particle displacement just after the plateau, a value 
that does not depend much on temperature. One possi- 
ble caveat is that this choice could yield poor "contrast" , 
i.e., only little differentiation between systems at differ- 
ent temperatures (we expect an increase in the patch 
correlation length as temperature decreases). 

To check the "contrast", consider two typical snap- 
shots of the system, one at low temperature and another 
at higher temperature, and evaluate the ratio of patch 
correlation lengths 7(0?*) = Ciowr/ChighT as a function 
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FIG. 1: Patch entropy as a function of patch volume: temper- 
atures range from T/T* = 4.35 (left) to T/T* = 0.59 (right) 
for a system of 2138 particles. V a is the volume of a single 
particle and T* the temperature of the avoided crystalliza- 
tion (see Ref. [13]). Note the three regimes of the entropy: 
constant and near zero for small sizes, then increasing lin- 
early, and finally saturation. The patch correlation length is 
extracted from the inverse of the slope of the linear (middle) 
regime, measured by a linear fit indicated by a dot-dash line 
for the lowest temperature. Insert: ratio of extracted length 
scales for T/T* = 0.35 and T/T* = 3.93 for the less fragile 
system versus the distance threshold involved in patch com- 
parison. In this article, the optimum value d* = 0.425 given 
by the maximum has been used. 



of d*. For d* — > 0, j(d*) -» 1, since in this limit of 
perfect resolution, each patch appears only once due to 
thermal fluctuations. On the other extreme, for d* — > oo, 
j(d*) — >• 1 as well, since in this limit all patches are con- 
gruent. For other values of d* , j(d*) > and is expected 
to have a single maximum - this is where the "contrast" 
is greatest. As seen in the inset to Fig. [T] a maximum is 
found at d* ~ 0.425cr (where a is the particle diameter); 
this value does indeed correspond to a mean displace- 
ment at a time just after the plateau [14] . The results we 
present here use this value of d* for a density per 2 ~ 0.85. 
In addition, we use an angular resolution 5a = 
which appears sufficient to approximate the 5a — s- lim- 
iting case well. 

Fig. [I] shows the patch entropy vs. patch size for var- 
ious temperatures. The entropy exhibits three regimes 
with increasing R: flat (subextensive) at small R, fol- 
lowed by a linear rise (the extensive regime), and finally 
saturation. The flat regime corresponds to small repet- 
itive patches, the second regime reflects the extensive 
behavior characteristic of a material disordered on this 
scale, and the saturation is a finite size effect (Increasing 
the system size moves the plateau higher). As discussed 
above, we extract the patch correlation length £s from 
the correlation volume given by the (inverse) slope of the 
extensive entropy regime, which, to average out oscilla- 
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FIG. 2: The structural correlation lengths vs. temperature 
for different fragilities: ^ (dark blue open symbols) denotes 
the hexatic length, £pts (light green open symbols) is the 
point-to-set length, and £s (solid red symbols) is the patch 
correlation length. Errors bars are roughly of the same size 
as the symbols. The vertical axis corresponds to absolute 
values of ^g. £pts and £s have been linearly scaled but their 
raw values are of the same order of magnitude as £pts is 
multiplied by a factor between 0.56 and 1.1 (depending on the 
curvature), £s by a factor between 0.52 and 0.65, and both 
have been vertically shifted. The three different curvatures k 
correspond to three different fragilities of the system, with the 
smallest curvature being associated with the largest fragility. 



tions due to the discrete particle size, we fit with a linear 
function (see Fig. [T]). We note that the behavior of the 
entropy for different temperatures clearly indicates that 
the order is growing as the system is cooled. 

In order to determine the relevance of the patch cor- 
relation length, it is of interest to compare it to other 
proposed measures. To this end, we have computed the 
point-to-set length £ PT s [H [71 [S] by adapting the method 
described in [3] to two-dimensional curved space [TH] • De- 
tails will be published elsewhere, but we find a slightly 
compressed exponential dependence of the overlap on the 
size of the cavity that is compatible with the results 
of [1]. Additionally, we have computed the hexatic or- 
dering length ^6 [11J ■ which serves as a benchmark for 
this system. 

We present the results in Fig. [2j where we plot the 
various lengths for different temperatures and fragili- 
ties. The three lengths have been linearly rescaled (see 
the figure caption). As can be seen, the three lengths 
show very similar temperature dependence for each of 
the three fragilities studied, except possibly at the lowest 
temperatures. The fragility dependence of these length 
scales (growing faster for larger fragilities) confirms the 
link between the cooperative character of relaxation and 
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FIG. 3: Density profile of the mean patch of the largest class 
of patches: the diameter is equal to 8.4<r, T/T* = 0.52, and 
na = 0.05. This class contains 1591 repeating patches out of a 
total of 4279 patches in the system, so roughly one third of the 
system is used to construct this density profile. Red (dark) 
spots indicate the most probable positions of particles inside 
of patch of this diameter, indicating a pronounced hexagonal 
order out to this scale. 

fragility. It is worth noting that the patch entropy 
method is in this case the least demanding in terms of 
computer resources, which is why it is computed down 
to lower temperatures. 

An additional advantage of the patch entropy method 
is that it can be used to extract information about the 
way in which the system orders. The idea is that the local 
order should be reflected in the largest classes of congru- 
ent patches, and the associated "mean patches" should 
exhibit the symmetries characteristic of patches in this 
class. To do this, we construct mean patches from the 
various equivalence classes computed for a given patch 
size from a given sample. This is done by superimposing 
and adding all the patches in the class (each at its opti- 
mal angle). This creates a density profile which reflects 
the symmetry of the class. Clearly this is relevant only 
if the population of the class is non trivial. In Fig. [3j we 
reconstruct the mean patch of the largest class of patches 
for a diameter 2R smaller than £s in the deeply super- 
cooled regime. The density profile clearly shows modula- 
tions indicating that the most probable particle locations 
are clearly arranged on a hexagonal lattice. For larger 
patches, the mean patch begins to blur, indicating that 
such patches are larger than the structural correlation 
length. From this, one could compute an absolute value 
of the correlation length by looking at the crossover be- 
tween peaked and blurred mean patches. For instance in 



Fig. [3j the patch begins to blur for a diameter just above 
8.4(7 (the one of the displayed patch). Thus £s ~ 9 — 10 a 
in this case, which seems inconsistent with results from 
Fig. [2j at this frustration and temperature £s ~ 25a. 
This may be rationalized by noticing that £s is, as stated 
above, defined up to a multiplicative constant and is lin- 
early rescaled to match £q m Fig. [2j showing a compara- 
ble temperature dependence but implying nothing about 
the absolute value of the lengths. 

Although local hexagonal order was expected for our 
system, these results demonstrate the feasibility of the 
patch entropy method as a way of identifying local or- 
der. It should hold equally well for other amorphous 
systems where the nature of the order is not known a 
priori. We are currently using this method to elucidate 
the symmetries of three-dimensional colloidal systems. 
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